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ABSTRACT 

The nearby neutron star low-mass X-ray binary, Cen X-4, has been in a quiescent state since its last outburst 
in 1979. Typically, quiescent emission from these objects consists of thermal emission (presumably from the 
neutron star surface) with an additional hard power-law tail of unknown nature. Variability has been observed 
during quiescence in Cen X-4 on both timescales as short as hundreds of seconds and as long as years. However, 
the nature of this variability is still unknown. Early observations seemed to show it was all due to a variable 
hard X-ray tail. Here, we present new and archival observations that contradict this. The most recent Suzaku 
observation of Cen X-4 finds it in a historically low state, a factor of 4.4 fainter than the brightest quiescent 
observation. As the spectrum during the brightest observation was comprised of approximately 60% from 
the thermal component and 40% from the power-law component, such a large change cannot be explained by 
just power-law variability. Spectral fits with a variable thermal component fit the data well, while spectral fits 
allowing both the column density and the power-law to vary do not, leading to the conclusion that the thermal 
component must be variable. Interestingly, we also find that the thermal fraction remains consistent between 
all epochs, implying that the thermal and power-law fluxes vary by approximately the same amount. If the 
emitting area remains unchanged between observations, then the effective surface temperature must change. 
Alternatively, if the temperature remains constant, then the emitting area must change. The nature of this 
thermal variability is unclear, but may be explained by variable low-level accretion. 
Subject headings: stars; neutron — X-rays; binaries — X-rays; individual; Cen X-4 



1. INTRODUCTION 

Neutron star low-mass X-ray binaries (LMXBs) are often 
transient. During outbursts they accrete at a significant frac- 
tion of the Eddington luminosity, yet these outbursts tend to 
only last a few months in most cases. Thus, for the major- 
ity of the time, these sources are in a quiescent state where 
the X-ray luminosity is significantly fainter. The X-ray spec- 
tra of neutron star LMXBs during quiescence usually con- 
sists of two components - a soft the rmal component and a 
harder power-law component (see e.g. lCampana et al]ll998aL 
for a review). The soft component is usually interpreted as 
thermal emission from the neutron star surface. The neu- 
tron star is heated during outbursts as compression by the ac- 
creted material causes nuclear reactions to occur deep in the 
crust. This heat is t hen radiated thermally during quiescence 
(iBrown et alJll998h . The hard power-law component is less 
well understood and may be associated with residual accre- 
tion, or pulsar shock emission (e.g., [Ca mpana et al. 1998b; 
ICampana & Stellal l2000'; 'Men ou & McClintockll20 1 ) . 

Accurately measuring neutron star radii is vital for discrim- 
inating between the large range of possibilit ies for the dens e 
matter equation of state (e.g., iLattimer & PrakashI 12007b . 
There are a number of potential methods for constraining 
neutron st ar radii, for instance using X-ray bursts (e.g., 
lOzel et all 1 2009; Guveretal. 2010), quasi-periodic oscilla- 
tions (j Miller et al. 1998) or relativistic Fe K emission lines 
dCackett et al.. .2008a) . However, one of the most promis- 
ing methods for measuring a neutron star radius, R, uses 
the thermal emission from quiescent neutron stars because 
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as with any blackbody-like emission / oc R^/D^ (where / 
is the source flux, and D the distance). While many quies- 
cent neutron stars are too faint for this method to produce 
constraining radius measurements, modest constraints have 
been possible for several o bjects (e.g., Heinkeetal. 2003] 
120061; I Webb & Barret! |2007|) . Future X-ray telescopes, such 
as the International X-ray Observatory, will allow accurate 
radius measurements from many quiescent neutron stars in 
the Galaxy. 

Nevertheless, there may be a potential problem. Qui- 
escent emission from neutron stars has been seen to be 
variabledC aiiipana et al.1 [19971 iRutledse et al.1 l200H l2002ah 
Campana & Stellal 120031 ICampana et al. 2004; Cackett et al] 
2005). In most cases this variability comes from observing 
sources at different epochs, and can be explained away by 
variations in the power-law, for instance due to changes in 
residual accretion rate (e. g., iCacke tt et al. 2005). However, 
in the case of Cen X-4 (C ampana et al. 2004), variability was 
observed during an XMM-Newton observation, and the na- 
ture of the variability remains unclear. If it is due to changes 
in temperature of the thermal component, this poses prob- 
lems for neutron star radius measurements. A mechanism 
for any short-timescale temperature change is uncertain, espe- 
cially under the standard deep crustal heating picture. Thus, 
this variability has important implications for radius measure- 
ments using thermal emission from quiescent neutron stars. 

1.1. Variability in quiescence 

Several explanations have been used to describe the X- 
ra y variability of quiesce nt neutron star LMXBs. In Aql X- 
1 iRutledge et al.l (l2002ah found that the temperature of the 
thermal component varied between observations. Conversely, 
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an analysis of the same data by ICampana & Stellal ( l2003h 
came to an opposing conclusion, suggesting that the variabil- 
ity could be attributed to correlated changes in the column 
density and the slope of the power-law. Moreover, an analysis 
of two observations of the neutron star LMXB in the globular 
cluster NGC 6440 led Cackett et al. (2005) to suggest that it 
was the power-law component varying in that object. 

Probably the best studied quiescent neutron star is Cen X- 
4, and multiple X-ray missions have observed it in quies- 
cence. Cen X-4 is one of the nearest known quiescent neu- 
tron star LMXB {D = 1.2 ±0.3 kpc; Chevalier et al. 1989). 
X-ray variability of Cen X-4 in quiescen ce has been known 
about for 10 years (ICampana et al.|[r997h . It has been seen 
to be variable over ti mescales of 5 years, varying by 40% 
dRutledge et al.l 12001*), o ver a period of 3 days, varying by 
a factor of ~3 (Campana et al.|[l997h and most recently on 
timescales as short as 100 seconds with an rms variability of 
45% (Campana et al. 2004). The spectrum of Cen X-4 clearly 
has both a thermal component and a power-law (with an in- 
dex in the range 1-2) present, yet, studies so far have been 
inconclusive as to whether it is the temperature of the thermal 
component, the spectral index of the power-law or the strength 
of the power-law component that is variable. 

Rut ledge et al.l (1200 Ih find a 40% change in the quiescent 
luminosity of Cen X-4 between ASCA and Chandra obser- 
vations ^5 years apart, which they attribute to the power- 
law component. However, on shorter timescales it is un- 
clear whether this interpretation holds. Campana et al. (2004) 
performed a detailed analysis of the quiescent spectrum of 
Cen X-4, finding variability on 100 s timescales throughout 
a ~50 ks XMM-Newton observation. A color-color analysis 
was not conclusive as to the source of the variability, thus, 
they extracted three separate spectra from the observation de- 
pending on the count rate. Based on an analysis of the low, 
medium and bright count rate spectra they remained unable to 
conclusively determine which component led to the variabil- 
ity. 

In this paper, we present an analysis of new Suzaku data 
of Cen X-4, as well as, analyzing two archival Chandra and 
XMM-Newton observations of this source. From spectral fit- 
ting, we show that the thermal component in Cen X-4 must be 
variable. 

2. DATA REDUCTION 

In Table [1] we detail all the observations analyzed here. Be- 
low we discuss the data reduction for each specific telescope. 

2.1. ASCA 

We use the pipeline produced spectra for the ASCA observa- 
tion, obtained from HEASARC. The re is a second AS CA ob- 
servation of Cen X-4 taken in 1997 dAsai et al.lll998l) . How- 
ever, as noted by 'Rutledg e et al.l (1200 Ih this observation is of 
significantly lower S/N than the first ASCA observation and 
therefore adds little extra to our aim of understanding the vari- 
abihty of Cen X-4. 

2.2. Chandra data reduction 

We reprocessed the data following the standard Chandra 
analysis threads, using CIAO version 4.2. Both observa- 
tions used the ACIS-S instrument. The psextract tool was 
used to extract the source and background spectra. For both 
Chandra observations we used a circular source region with 
a 10 pixel radius, and for the background an annulus with 



inner radius of 15 pixels and outer radius of 45 pixels cen- 
tered on the source. The response files were created using 
the mkacisrmf and mkarf tools. The source spectra were 
binned to a minimum of 20 counts per bin in the energy range 
0.3 - 10 keV to allow the use of x^ statistics in spectral fitting. 

2.3. XMM-Newton data reduction 

Data were analyzed using the XMM-Newton Science Anal- 
ysis Software, version 9.0.0. Calibrated event lists were cre- 
ated from the Observation Data Files using the latest calibra- 
tion files. The first XMM-Newton observation was performed 
with the MOS and PN detectors operated in prime full win- 
dow mode with the thin filters. This first XMM-Newton obser- 
vation suffered large background flaring throughout the obser- 
vation. We therefore only extracted spectra from times with 
low background. For the MOS detectors, we excluded times 
when the count rate from events > 10 keV was higher than 0.5 
c/s. For the PN detector, we excluded times when the count 
rate from events in the range 10-12 keV was higher than 
1.0 c/s. This reduced the net exposure times to 36.7, 36.5 and 
23.0 ks for the MOS 1, MOS 2 and PN detectors, respectively. 
We used a circular source extraction region of radius 45" for 
MOS detectors and 40" for the PN. For background extrac- 
tion regions we used a source-free, nearby 2' circular region 
for the MOS detectors, and two source-free, nearby circular 
regions of 40" and 60" for the PN. Net count rates were 0.27, 
0.29 and 1.2 c/s for the MOS 1, MOS 2 and PN, respectively. 
All source spectra were binned to a minimum of 20 counts per 
bin. 

The second XMM-Newton observation was performed with 
the MOS and PN detectors operated in timing mode. This sec- 
ond XMM-Newton observation also suffered from significant 
background flaring, therefore we filtered the data, excluding 
times with high background. For the MOS detectors we ex- 
cluded times when the count rate from events > 10 keV was 
higher than 0.05 c/s. For the PN detector, we excluded times 
with the count rate from events in the range 10-12 keV was 
higher than 0.5 c/s. The net exposure times were reduced to 
64.3, 63.4 and 58.9 ks for the MOS 1, MOS 2 and PN, re- 
spectively. The source extraction regions had RAWX values 
of 306 - 330, 296 - 320 and 31 - 45 for the MOS 1, MOS 
2 and PN. Background extraction regions were taken towards 
the edge of the exposed part of the detectors, away from the 
source, and had RAWX values 260 - 284, 261 - 272 and 7 
- 21 for the MOS 1, MOS 2, and PN. Net count rates were 
0.19, 0.19 and 0.86 c/s (MOS 1, MOS 2 and PN). The back- 
ground rate was higher during this observation than the first, 
thus significantly higher counts per bin was required to give 
significant detections in all bins. We binned to a minimum of 
150 counts per bin for MOS 1 and MOS 2, and 500 counts per 
bin for the PN. 

2.4. Suzaku data reduction 

The Suzaku observations were performed with a full win- 
dow, and with normal clocking modes, with the telescope at 
the nominal XIS aimpoint. The data were analyzed using 
HEASoft version 6.8, which includes the Suzaku vl5 soft- 
ware. 

Data reduction for the XIS detectors, which cover the soft 
X-ray energy band (approximately 0.5 - 10 keV), follows. 
The source spectra were extracted for each detector (XIS 0, 
1 and 3) using a circular region of radius 250 px, centered 
on the source. The background spectra were extracted from 
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TABLE 1 
Observations of Cen X-4 



Sliort name 



Start date 



Exposure time(ks) 



Reference 



ASCA 

Chandra 

XMM-Newton 

XMM-Newton 

Chandra 

Suzaku 



41008100 

713 

0067750101 

0144900101 

4576 
403057010 



ASCA 
CXOl 
XMMl 
XMM2 
CX02 
SUZ 



27/02/1994 
23/06/2000 
20/08/2001 
01/03/2003 
21/06/2004 
16/01/2009 



39 
10 

53 
78 
10 

147 



1,2,3,4 
4 
5 
6 
6 
6 



References. — (D IAsai et al1[T99a (2) IAsaiet all 19981 OI IRutledgeet all 19991 (4) IRutledge et aU2001L (5) ICampana et aU2004 (6) This work 
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Fig. 1. — Unfolded spectra from the brightest (XMMl, black) and faintest 
(SUZ, red) quiescent observations of Cen X-4. The neutron star atmosphere 
(blue, dashed line) and power-law (green, dotted line) components are also 
shown. The large variability requires that the thermal component has to vary 
between the two epochs. For XMMl we only show the PN spectrum, and for 
SUZ we show just the combined XIS -H 3 spectrum. The spectra have been 
rebinned for visual clarity. 



an annulus with inner radius 300 px and outer radius 425 px, 
centered on the source. The responses were generated with 
the xisrmfgen and xissimarfgen tools. We co-added 
the spectra from the two front-illuminated detectors (XIS 
and XIS 3) to increase the signal-to-noise ratio. The spectra 
we binned to a minimum of 250 counts per bin in the 0.5 - 10 
keV energy range. 

We also extracted a spectrum from the hard X-ray detec- 
tor PIN camera, which covers the energy range from approx- 
imately 10 - 70 keV. The spectrum was extracted using the 
hxdpinxbpi tool which also extracts the background spec- 
trum. The net source exposure was 124 ks, however, the 
source was not detected, and the spectrum is consistent with 
the background. 

3. SPECTRAL ANALYSIS 

The spectrum of Cen X-4 has clearly been shown in the 
past to have both a thermal and non-thermal component 
(see Fig. [T]i- We therefore choose to fit the spectra with 
an absorbed neutron star hydrogen atmosphere plus power- 
law model. We use the nsatmos model (iHeinke et alJ 
12006) which includes thermal electron conduction and self- 
irradiation by photons from the compact object. It also as- 
sumes a negligible magnetic field (< lO** G), which is relevant 
here. Spectral fitting is performed using XSpec vl2 (lArnaudI 
119961) throughout, and all uncertainties quoted are at the Icr 
level of confidence. 



The aim is to investigate the variability of Cen X-4. 
Through spectral fitting we test whether the thermal compo- 
nent is variable, or whether changes in the absorbing column 
and power-law component can explain the variability. We 
therefore fit the spectra simultaneously, tying several param- 
eters between the spectra and letting others vary freely. We 
do this in several ways. Firstly, we allow the temperature of 
the neutron star atmosphere to vary between observations, but 
hold the emitting radius fixed at 10 km. Secondly, we hold 
the temperature tied between the observations (i.e. the tem- 
perature is the same for all observations, but this temperature 
is a free parameter in the fit) but allow the emitting radius to 
vary. In both cases the power-law component is left free to 
vary, and the photoelectric absorption column density (XSpec 
model phabs) is tied between all observations. Finally, we 
also test the possibility that there are changes in both the col- 
umn density and the power-law while the thermal component 
remains unchanged. In this case both the absorption column 
density, A^h and the power-law component (both the index 
and normalization) are free to vary between epochs, while the 
temperature and radius of the neutron star atmosphere com- 
ponent are tied between the epochs (but still a free parameter 
in the fit). All spectral fits are within the energy range 0.5 - 
lOkeV. 

The ASCA observations consist of four spectra (2 GIS and 
2 SIS). The parameters are tied between all the A5CA spectra. 
Similarly for the Suzaku observation we have two spectra, one 
from the combined XIS H- 3, and one from the XIS 1 detec- 
tor Again, parameters are tied between the detectors. For 
the first XMM-Newton observation we have three spectra (two 
MOS, one PN), and tie all parameters between the detectors. 
When fitting the second XMM-Newton observation on its own, 
we found a offset between the MOS and PN detectors (which 
were both operated in timing mode). Fitting the absorbed neu- 
tron star atmosphere plus power-law model to all three spectra 
with all parameters tied between the detectors led to a poor fit 
(X^ =2.81). However, if a constant offset was used an accept- 
able fit was achieved (x^ = 1.02). The constant was fixed to 
1.0 for the PN, and the best-fitting value for the MOS spectra 
was 1.26 ±0.01 (the constant was tied between the MOS 1 
and MOS 2). Comparing fitting the PN alone and fitting all 3 
spectra with a constant, we found only very minor differences 
in the fit parameters. Thus, in the remaining analysis, we fit 
all 3 spectra with a constant offset between the PN and MOS 
to achieve the best constraints on parameters. 

Spectral fit results from allowing the temperature to vary are 
given in Table |2l spectral fit results from allowing the radius 
to vary are given in Table [3] and spectral fits from allowing 
both the column density and power-law to vary are given in 
Table S 
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TABLE 2 
Spectral fits for Cen X-4, with temperature variable and radius fixed 



ASCA 



CXOl 



XMMl 



XMM2 



CX02 



suz 



kT^^ (eV) 
Power-law index, 
Power-law norm. 



r 

(10-- 



cm-2) 



Unabs. 0.5 
Xi (dof) 



10 keV thermal flux (10"'^ erg s"' cm^^) 



63.3 ±0.7 


59.2 ±0.4 


4.9: 
66.2 ±0.2 


±0.2 

62.0 ±0.2 


51.2±0.6 


48.2 ±0.6 


1.24 ±0.17 


0.97 ±0.25 


1.41 ±0.05 


1.26 ±0.08 


0.78 ±0.43 


1.69±0.17 


7.9±1.6 


3.1 ±0.9 


9.4 ±0.6 


5.4 ±0.5 


1.4 ±0.7 


3.7±0.6 


1.86±0.17 


1.20±0.13 


2.09 ±0.06 


1.48 ±0.03 


0.63±0.10 


0.47 ±0.03 


1.00 ±0.09 


0.71 ±0.07 


1.26 ±0.04 


0.90 ±0.02 


0.33 ±0.06 


0.24 ±0.02 






0.99(1166) 







Note. — The neutron stai' atmosphere model 'nsatmos' is used. For all observations, the neutron star radius was fixed at 10 km, and the mass at 1.4 Mq. The 
distance was fixed at 1 .2 kpc. Wh was tied between all observations. kT^ is the effective temperature for an observer at infinity, for i? = 10 km and M = 1 .4Mq . 
The power-law normalization is defined as photons keV"' cm"^ s"' at 1 keV. 



TABLE 3 
Spectral fits for Cen X-4, with temperature tied and radius variable 





ASCA 


CXOl 


XMMl 


XMM2 


CX02 


SUZ 


Wh (lO^o cm-2) 






5.6±0.6 






kT,ff (eV) 






79.3 ±2.5 






Radius (km) 


10.8 ±0.6 


9.7 ±0.6 


11.7±0.8 


10.5 ±0.5 


7.8 ±0.4 


7.2±0.3 


Power-law index, F 


1.29±0.11 


0.85 ±0.14 


1.51 ±0.04 


1.22 ±0.05 


0.31±0.16 


1.53±0.10 


Power-law norm. (10^^) 


8.5 ±1.7 


2.6 ±0.8 


11.1±0.8 


5.1 ±0.5 


0.8 ±0.5 


3.0 ±0.6 


Unabs. 0.5 - 10 keV flux (10-'^ erg s"' cm-2) 


1.88±0.15 


1.24±0.12 


2.13±0.05 


1.51 ±0.05 


0.72±0.10 


0.48 ±0.04 


Unabs. 0.5 - 10 keV thermal flux (lO^'^ erg s"' cm-^) 


1.01 ±0.08 


0.74 ± 0.07 


1.27 ±0.03 


0.94 ±0.03 


0.35 ±0.05 


0.25 ±0.02 


xi (dof) 






0.99(1165) 







Note. — The neutron star atmosphere model 'nsatmos' is used. For all observations, the mass was fixed at 1.4 M0. The distance was fixed at 1.2 kpc. Both 
Nn and kT^ff were tied between all observations. Note that here we give the unredshifted kT^g as the emitting radius (and thus the redshift) is allowed to be 
different at each epoch. The power-law normalization is defined as photons keV"' cm"- s"' at 1 keV. 



TABLE 4 
Spectral fits for Cen X-4, with variable A'h and power-law 





ASCA 


CXOl 


XMMl XMM2 


CX02 


SUZ 


Nh (lO^o cm-2) 

kT^^ (eV) 

Radius (km) 

Power-law index, F 

Power-law norm. (10^^) 

Unabs. 0.5 - 10 keV flux (lO^'^ erg s"' cm-^) 

xi (dof) 


12.0 ±1.4 

2.0±0.I 

20.1 ±2.2 
2.10±0.09 


14.6±0.8 

1.3±0.1 

5.1±1.0 

1.61±0.11 


10.0±0.6 11.6±0.7 

48 7+''-5 

25.7 ±0.5 

2.4±0.1 1.9±0.1 

32.1 ±1.2 13.1 ±0.9 

2.31 ±0.04 1.78 ±0.03 

1.47(1165) 


31.5±1.5 

-1.0+«'^ 
0.16+"'^ 
1.99+?|] 


44.9 ±1.7 

0.70±0.14 
0.87 ±0.23 

, ^9+0. 14 



Note. — The neutron star atmosphere model 'nsatmos' is used. For all observations, the mass was fixed at 1.4 M0. The distance was fixed at 1.2 kpc. Both 
the column density, N^, and the power-law (index and normalization) were allowed to vary between epochs, while the neutron star effective temperature and 
radius were tied between all observations. kT^ is the effective temperature for an observer at infinity, for M = 1 .4Mq and the best-fitting radius. The power-law 
normahzation is defined as photons keV"' cm"^ s"' at 1 keV. 



4. LIGHTCURVES 

In addition to studying the long-term variability of Cen X- 
4 from epoch to epoch, we also take a look at the variability 
during the observations newly analyzed here (XMM2, CX02, 
SUZ). We extracted background-subtracted lightcurves using 
the same source and background regions as for the spectral 
analysis. Given the low count rate, the lightcurve from the 
Suzaku observation is noisy, and does not show any clear 
variability. However, the lightcurves from both XMM2 and 
CX02 show some short-term variability. The lightcurves 
from XMM2 are shown in Fig. |2] and the lightcurve from 
CX02 is shown in Fig. [3] 

The XMM2 lightcurves were extracted from the entire data, 
including the time periods that were excluded from the spec- 
tral analysis due to high background flaring. We show the 



background-subtracted 0.3 - 10 keV lightcurves with 250-s 
binning from all three EPIC detectors, as well as showing the 
background lightcurve from the PN. Some of the variability 
seen in the lightcurve may be associated with imperfect back- 
ground subtraction during periods of high background flaring. 
A clear example of this is at around 20 ksec. However, there 
is a significant flare in the source lightcurves at around 37 
ksec that is during a period of low background, thus is clearly 
associated with source variability. 

In order to investigate the nature of this flare, we extracted 
lightcurves in the 0.3 - 2.0 keV range and the 2.0 - 10 keV 
range from the PN data, and looked at the hardness ratio (2.0 
- 10 keV count rate/ 0.3 - 2.0 keV count rate) for any sig- 
nificant changes (see Fig. |4|. The flare is prominent in the 
soft lightcurve, though there is no significant evolution in the 
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Fig. 2. — Lightcurves from the second XMM-Newton observation (XMM2). The top three panels show the MOS 1, MOS 2 and PN net count rate, with 250 
s binning in the 0.3 - 10 keV energy range. The bottom panel shows the PN background count rate (note that periods of high background excluded from the 
spectral analysis have not been filtered out here). While some variability in the net lightcurves may be associated with higher levels of background (for instance 
at around 20 ksec), there is a clear flare from the source at around 37 ksec seen in all three detectors during a period of low background. 
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Fig. 3. — 0.3 - 10 keV lightcurve from the second Chandra observation 
(CX02), with 250 s binning. The dashed line is the weighted average of the 
lightcurve. The early section of the lightcurve is significantly higher than the 
average. 



hardness ratio during the flare. It is therefore not possible to 
conclude which component caused the flare. 

We perform a simple test for variability in the CX02 obser- 
vation by fitting a constant to the lightcurve. A constant pro- 
vides a poor fit (x^ = 3 .0), with the early part of the lightcurve 
significantly above the weighted average (see Fig.O, indicat- 
ing variability on a timescale of a few hundred seconds. 

5. RESULTS 

The spectral analysis clearly shows a large amplitude of 
variability between epochs, with the flux varying by a factor 
of 4.4 between the brightest and faintest observations (XMMl 
and SUZ, respectively). We show the unfolded spectra from 
these two observations in Figure[Tl and the spectral fits shown 
are for varying temperature and fixed radius. What is clear 
from the spectral fits is that the thermal component has to 
have varied between these two epochs. Given that the power- 
law accounts for approximately 40% of the 0.5 - 10 keV flux 
in XMMl, a factor of 4.4 change in flux cannot be achieved 
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Fig. 4. — The XMM2, PN lightcurves (with 250 s binning) around the time 
of a significant flare. Here we show the 0.3 - 2 keV (top) and 2-10 keV 
(middle) lightcurves, and the associated hardness ratio (2-10 keV count rate 
/ 0.3 - 2 keV count rate). There is no obvious evolution in the hardness ratio 
during the flare. 
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through a variable power-law alone. The spectral fits allowing 
both the power-law and column density to vary also do not fit 
the data well, leading to a significantly worse fit (x^ = 1 .47, 
V = 1165) than when allowing the thermal component to be 
variable. 

A spectral fit with the temperature tied between all obser- 
vations, and the radius fixed at 10 km is statistically not ac- 
ceptable (x^ = 1.91, v=\ 171). The brightest observations all 
have very soft power-law indices (^ 3), and poor fits above 3 
keV, as the power-law component tries to fit the majority of 
the thermal component. We also tried fitting with the radius 
as a free parameter, though having its value be the same for all 
observations, however this does not significantly improve the 
fit (x^ = 1.90, V = 1170). Thus, the spectral fits demonstrate 
that the thermal component must vary between epochs, and 
the variability cannot be attributed to changes in the power- 
law index, normalization and/or the absorption. 

Figure |5] shows the long-term variability of Cen X-4. The 
top three panels in this figure show the variability of the (a) 
unabsorbed 0.5 - 10 keV flux, (b) effective temperature and 
(c) fraction of flux in the thermal component from the spectral 
fits where the temperature is variable between epochs, and the 
radius fixed. Both the changes in flux and temperature follow 
the same overall pattern. In fact, it can be seen that the thermal 
fraction remains consistent between epochs, implying that the 
flux from the thermal and power-law components is varying 
by approximately the same amount. The bottom panel, (d), 
shows how the radius changes in the case where the tempera- 
ture is tied between each epoch in the spectral fits. It demon- 
strates that in this model, significant changes in the radius are 
required to account for the spectral changes. Therefore, to 
achieve the observed flux variability, either the temperature 
of the thermal component, or the emitting radius must vary 
between epochs. We cannot rule out that both the temperature 
and the radius change. 

For the sake of completeness, we also note several pre- 
vious quiescent observations not studied here. Einstein and 
EXOSAT o bservations were perform ed in 1980 and 1986 re- 
spectively (Ivan Paradijs et al.lll987) . though their flux mea- 
surements ar e not precise enough t o be constraining here (see 
table 2 from iRu tied ge et al]|200lL which concisely summa- 
rizes these observations). R OSATfHRl observed Cen X-4 in 
1995 dCampana et al.|[T997l) . however, given that this was not 
a spectral instrument no tight flux constraint is obtained. As 
noted earlier, the re was a second A SCA observation of Cen X- 
4 in quiescence (lAsai et alJI 19981) . which was unfortunately 
of low quality. This observation was perfo r med o n Febru- 
ary 4-5, 1997. The analysis by lAsai et all (119981) led to a 
rather poorly constrained unabsorbed 0.5 - 10 keV flux of 



;+1.75 
-1.17 



1.75;!;r,T x lO erg s cm , consistent with all but our 



faintest observation. Finally, a BeppoSAX observatio n was 
performed on February 9, 1999 (ICampana et alJl200Q') . and 
this wa s found to be consist ent with the first Chandra obser- 
vation ( lRutledgeetalJl2001h . 
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Fig. 5. — (a) Long-term lightcurve of Cen X-4. The data shows the unab- 
sorbed 0.5 - 10 keV flux from fitting with the temperatures tied and a radius = 
10 km. (b) Effective temperature (for an observer at infinity) vs. time. In the 
spectral fits the neutron star radius was fixed at 10 km. (c) Thermal fraction 
vs. time, where the thermal fraction is defined as the ratio of the unabsorbed 
0.5 - 10 keV neutron star atmosphere flux to the total unabsorbed 0.5 - 10 
keV flux. This is using the spectra fits with the temperature tied and a radius 
= 10 km. (d) Neutron star radius vs. time from spectral fits assuming that 
the neutron star atmosphere temperature is the same at all epochs and that the 
radius varies. 



6. DISCUSSION 

The two known outbursts from the neutron star LMXB 
Cen X-4 were in 1969 and 1979, and the source has been 
in quiescence ever since (it has not been detected in outburst 
by any mission or all-sky monitor). Our spectral analysis of 
six observations during quiescence covers a period of over 15 
years, and shows variability of a factor 4.4 between the bright- 
est and faintest states. In all observations both a thermal and 



power-law component is present in the spectrum. The am- 
plitude of the observed variability requires that the power- 
law alone cannot account for the variability, as the power- 
law component only contributes about 40% of the flux dur- 
ing the brightest state. Our spectral fitting demonstrates that 
the thermal component must vary with either the temperature 
and/or the radius of the thermal component changing between 
epochs. Allowing both the power- law component and the col- 
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umn density to vary (while the thermal component remains 
constant) does not fit the data well. 

In addition to the long-term epoch to epoch variability, we 
have observed short-term variability on timescales of the or- 
der of a few hundred seconds in the lightcurves from two 
of the newly analyzed datasets (XMM2 and CX02). Such 
sh ort-term variabihty ha s previously been observed in Cen X- 
4 JCampana et al.ll2004l) . though the nature (whether thermal 
or power-law) was inconclusive. Given the source was fainter 
during these two new observations than during the study by 
ICampana et al. (2004), we again are unable to determine the 
cause of this short-term variability. 

There are several methods by which the temperature of the 
neutron star surface may change on timescales of years, which 
we now discuss in turn. Firstly, once a source returns to qui- 
escence, the crust may cool (e.g. Rutledge et al. 2002b). In 
the case where the outburst is particularly long (i.e., lasting 
several years rather than several weeks to months), the crust 
can be heated significantly out of thermal equilibrium with the 
rest of the star. Therefore, once accretion reduces to quiescent 
levels, the crust will thermally relax. Such crustal cooling has 
been observed in four sources so far. It was first observed 
in KS 1731-260 and MXB 1659-29 which both returned 
to quiescence in 2001 dWiinands et all 1200 IL I2003h . Since 
then, monitoring with Chandra and XMM-Newton has shown 
both sources to co ol rapidly dWiinands et alj 120021 12004 
ICackettetal]|2006h . The cooHng curve of MXB 1659-29 
covers 6.6 years, and is well described by an exponential de- 
cay to a constant level, with a e-folding timescale of 465 ± 25 
days ( Cackett et al. 2006, 2008 b). While the cooling curve of 
KS 173 1-260 can also be described by e xponential cooUng t o 
a constant, with a similar e-folding time dCackett et alj|2006l) . 
the most recent observation suggests that cooling is continu- 
ing following a power-law decay (Cackett et al., in prepara- 
tion). 

In the last couple of years, two addition al sources with 
long outbur sts have also gone i nto quiescence dPegenaar et alj 
|2009; Frid riksson et aljT2010[) . The neutron star transient 
EXO 0748-676 was in outburst before returning to quies- 
cence. Chandra and Swift observations covering the first 5 
months after the source transitioned to quiescence shows an 
initially slow decrease in temperature (D egenaar et al. 2009 ), 
with further observations ongoing. Finally, once the transient 
XTE J1 701-462 returned to q uiescence it displayed rapid 
cooling dFridriksson et al.|[2010l) . However, the cooling curve 
observed is complicated by a temporary increase in tempera- 
ture about 220 days into quiescence. After this increase the 
cooling continued on the same track as before the increase 
- an exponential de cay to a constant, with an e-folding time 
of about 120 days. iFridriksson et alj d2010l) suggest that the 
apparent increase in temperature could be caused by an in- 
creased level of accretion. 

The long-term quiescent lightcurve for Cen X-4 (Fig. |5]l 
may show an overall decrease suggestive of cooling, but with 
several increased points, as observed in XTE J1701-462. 
Clearly, though, the sampling rate is extremely infrequent 
and there are only a small number of points. However, the 
timescale of any apparent decrease is significantly longer than 
the e-folding timescales for the crustal cooling sources dis- 
cussed above. The variability we observe from Cen X-4 is 
from 15-30 yr after the 1979 outburst, and thus the crust 
should have cooled back into thermal equilibrium with the 
core at this point if it has a similar structure to the other 
sources. It therefore seems unlikely that crustal cooling con- 



tributes significantly to the variability observed. Note that the 
neutron star core wi ll not cool appreciably over the timescales 
observed here (e.g..! Yakovlev & Pethickll2004l) and thus will 
not contribute to any variabihty. 

Continued low-level accretion onto the neutron star sur- 
face can also change the thermal properties of the star. 
The thermal quiescent flux is sensitive to both the amount 
of H/He remaining on the surface after an outburst and 
on the composition of ashe s from previous H/"*He burning 
dBrown. Bildsten. & Chan g "2002). Continued low-level ac- 
cretion will change the surface composition by adding an in- 
sulating layer of ashes from H/"*He burning. This therefore 
changes the surface effective temperature, and can occur on 
timescales of > 10 yr at the luminosity of Cen X-4. However, 
this effect should lead to an increase in temperature over time. 
Specific calculations for Cen X-4 by Brown et al. (2002) sug- 
gest only a 20% increase in brightness over 30 years due to 
this effect. Here, we see variations greater than 20% and vari- 
ability that is not just a simple increase, suggesting that an 
increase in the depth of the H/"^He layer cannot explain the 
observed variability in Cen X-4. 

Another process by which the thermal properties of a qui- 
escent neutron star may vary is diffusive nuclear burning 
dChang & Bildstenll2003h . In this process, hydrogen diffuses 
to deeper layers in the crust where it can fuse. This would 
lead to a drop in flux as hydrogen is consumed and the hydro- 
gen abundance in the photosphere decreases. This process is 
only exp ected to lead to a changes on timescales > 10 yr For 
Cen X-4. IChang & Bildstenl d2003h specifically predict that 20 
yr after the outburst the flux will vary by only 3% over a 10 yr 
timescale, and that roughly 100 yr after the outburst the flux 
will drop by 12% over a 20 yr period. Consequently, it is hard 
to explain the observed variability by this process, given that 
we see variability of much greater amplitude than this, as well 
as more than just a simple flux decline. 

The above explanations are also disfavored because the 
thermal fraction is constant while the luminosity varies. How- 
ever, one other possibility is that variable residual accretion is 
responsible for the observed behavior. If low-level accretion 
onto the neutron star surface is causing the power-law com- 
ponent in the spectrum, then an increase in this accretion rate 
would increase the power-law, which could in turn lead to an 
increase in the observed neutron star temperature. In fact, 
residual accretion is known to be able to produce a thermal 
spectral shape (Zampieri et al. 1995). One striking finding is 
that while the total luminosity varies by a factor of over 4 
the ratio of thermal to power-law flux remains approximately 
constant (Fig. |5] panel c). This appears to indicate that the 
thermal and power-law components are linked together, thus 
any mode for residual accretion would need to produce both 
components and in a way that their ratio remains constant. 
This may be able to provide constraints on models of accre- 
tion flows at low mass-transfer rates. 

In summary, the crustal cooling, H/He ashes and diffusive 
nuclear burning processes all seem unlikely to drive the ther- 
mal variability observed, and the most likely scenario appears 
to be one where variable low-level accretion causes both the 
thermal and non-thermal variability. More frequent monitor- 
ing of Cen X-4 will lead to a better understanding of both 
the amplitude and timescales of the variability. Finally, this 
thermal variability may pose a problem for measuring neu- 
tron star radii from modeling quiescent emission given that it 
is not clear whether the emitting radius, temperature or both 
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is causing the variability. If the variabiHty is due to low-level 
accretion, then avoiding sources with a significant power-law 
component would be beneficial. 
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